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Abstract 

We calculate the equation of state of nuclear matter in the self-consistent T-matrix 
scheme including three-body nuclear interactions. We study the effect of the three- 
body force on the self-energies and spectral functions of nucleons in medium. 
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The extrapolation of the energy per particle in dense nuclear systems from 
finite nuclei to infinite nuclear matter and to neutron matter in neutron stars 
involves theoretical estimates. Most of the approaches are using free nucleon- 
nucleon potentials when treating systems with many particles. Obviously a 
number of approximations must be done to obtain a result. Properties of 
dense nuclear matter have been estimated with a variety of different schemes : 
Brueckner-Hartree-Fock (BHF) calculations [Tf2"f3"f4"] . variational calculations 
[5l6lTTll8ll9] and self-consistent T-matrix calculations [T0|lll|ll2|ll3|ll4j . The cen- 
tral issue of these studies is the calculation of the binding energy of cold 
symmetric and neutron rich nuclear matter. It has been realized that a realis- 
tic description of the nuclear matter at saturation density and beyond cannot 
neglect the presence of three-body forces between nucleons. Variational and 
BHF calculations that take the three-body forces into account reproduce the 
empirical saturation point density and the binding energy in symmetric nu- 
clear matter p3|16lll7|l5|6lfT8] . The spectral T-matrix method using in medium 
dressed nucleon propagators and the self-consistent T-matrix is well suited for 
the calculation of single-particle properties of nucleons in the medium [TO"] , 
the effective scattering [19] and pairing correlations [2"0~|2~T] . One important 
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feature of the scheme is the automatic fulfillment of thermodynamic consis- 
tency relations |22j . The equation of state of symmetric and neutron matter 
in the T-matrix approach have been evaluated as well [12j2 5fTTl[T5] . however 
only two-body interactions have been included up to now. The resulting equa- 
tion of state is similar to the one obtained by other methods with the same 
interactions. In this letter we present first results for the self-consistent T- 
matrix approximation with phenomenological three-body interactions taken 
into account. After reproducing the properties of the symmetric nuclear mat- 
ter around the empirical saturation point we calculate the modified in-medium 
properties of dressed nucleons at several densities. 

The in-medium self-consistent T-matrix [2~if2"o1 26] is defined by the two-body 
potential V as 

T = V + VGGT (1) 



which denotes the resummation of the two-nucleon ladder diagrams. The in- 
medium nucleon propagator 

° = u-p 2 /2m-Y, ^ 



is dressed by the self-consistently calculated self-energy 

iT, = Tr [TG] . (3) 



The numerical calculations are performed in the real-time formalism for the 
finite-temperature Green's functions [10]. The spectral function for dressed nu- 
cleons is obtained from the retarded self-energy (with energy and momentum 
arguments explicitly written) 

A{p u) — —2\vaL(jp,uj) 

{uj — p 2 /2m — Re£(p, a;)) 2 + ImS(p, oj) 2 



At each density p the above set of equations is iterated until convergence and 
the Fermi energy /i is adjusted to fulfill the constraint 



du f d 3 p 



2tt J (2tt] 



A(p,u). (5) 



The binding energy can be calculated from the Galitskii-Koltun's sum rule 

E 1 r duj f d 3 p 

1 1 (LO + p 2 /2m)A(p,uj) , (6) 



N 2pJ27vJ (2tt 
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Fig. 1. The equation of state for symmetric nuclear matter. The solid line represents 
the result of the T-matrix calculation, the dashed line the BHF result [TTJ and the 
dashed line the variational result [6j. The lines with symbols denote the results of 
calculations including only the two-body force. 

a formula that works for Hamiltonians with two-body interactions. In the gen- 
eral case one should calculate the diagrams corresponding to the expectation 
value of the Hamiltonian < H > [27]. 

Several parameterizations of the nuclear two and three-body interactions are 
used in nuclear matter and finite nuclei calculations. Since the short range be- 
havior of the nuclear force is not precisely known, differences in parameteriza- 
tion of the nucleon-nucleon interaction at high momenta can be compensated 
by differences in the associated three-body term [28.29J. Explicit calculations 
using relativistic BHF formalism require a different (if any) three-body force. 
Still other three-body interactions are needed when using renormalized effec- 
tive two-body potentials [301 • ^ n the following we use a simple phenomeno- 
logical way of taking the three-body interaction into account [31] and for the 
two-body potential we take the CD-Bonn interaction. The three body term 
motivated by the two-pion exchange process has the form of an additional 
density dependent two-body interaction 
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V(r) = I c T 2 (r) (e-^ p - l) (7) 



where 



T(r) 



( 
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f] = 0.7fm _1 , c = 2fm™ 2 , I c = — 5.7MeV, 71 = 0.15fm 3 . The density dependent 
part is short range and repulsive and leads to a stiffening of the equation of 
state. We perform iterative self-consistent calculation of the dressed propaga- 
tors and the in-medium T matrix with such density dependent interactions at 
several densities between p = 0.6po and p = Spo, where the saturation den- 
sity is pQ = 0.17fm -3 . Additionally a phenomenological attractive mean field 
energy is taken in the form [31] 



where 72 = — 260MeV fm 6 and 73 = llfm 3 . With the density dependent two- 
body interaction (JTj) the Galitskii-Koltun's sum rule gives the same binding 
energy as the expectation value of the Hamiltonian. For the chosen interaction 
the difference shows up only in the term Etna, that should be added to the 
energy but the corresponding shift in the Fermi energy is 5(pEtna{p)/N)/5p. 

The binding energy in symmetric nuclear matter is shown in Fig. [I] and com- 
pared to results of variational and BHF calculations with two-body and three- 
body forces [T7]l6] . We observe that the inclusion of density dependent forces 
suffices to reproduce the empirical saturation point of nuclear matter, although 
our equation of state is slightly stiffer with the parameters of the short range 
density dependent interaction taken directly from |31j . As expected, in order 
to get realistic binding energies at higher densities three-body forces must be 
included in the calculation. The equation of state fixes the strength of the 
additional three body terms. 

Within the T-matrix scheme we explore the consequences of the additional in- 
teraction terms on the single-particle properties. In Fig.[2]is shown the spectral 
function A(p, u) at po- We see that the modification of the force changes the 
spectral function. This indicates that besides the equation of state the spectral 
function is sensitive to the chosen nuclear interaction. The binding energy is 
not sensitive enough to constraint the short range part of the nuclear inter- 
action. The nuclear spectral function when compared to experimental results 
[32] gives insight into the nuclear two and three-body interactions in nuclei. In 
a future publication we shall present the results of an investigation of the sen- 
sitivity of the proton spectral functions on the density and isospin dependence 
of the three-body force. Within the parametrization of the three-body terms 



E TNA (p)/N = 3j 2 p 2 e 



(9) 



4 




Fig. 2. The spectral function for dressed nucleons at saturation density. The dashed 
line is the result with the three-body force, the solid line represents the result of the 
calculation the with two-body interactions only. 

taken in this letter, the spectral functions shows a stronger quasi-particle peak, 
which is a consequence of a reduced scattering at the Fermi surface. 

It is interesting to analyze the real-part of the self-energy ReS(p, uj p ) at the 
quasi-particle pole (u p = p 2 /2m + ReS(p, ui p )). It is the in medium potential 
felt by the nucleon. The self-energy is given by dispersion relation 



ReE(p,uj) = E HF (p) + J 



duo' Im£(j>, 10 ) 



7T UJ — UJ 

^hf{p) + S disp (p,u;) 



(10) 
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Fig. 3. The real part of the self-energy at the quasiparticle pole Re£(p, u p ), at 
the saturation density po (upper panel) and at 3po (lower panel). We show the 
Hartree-Fock contribution, the dispersive part, and the sum (Eq. [10]) . The dashed 
lines are the results including three-body forces and the solid lines represent the 
energies obtained with two-body forces only. 

It is a sum of a dispersive self-energy £<&sp and the mean-field self-energy 
Y,hf(p). Different parameterizations of the nuclear interactions yield different 
Hartree-Fock and dispersive self-energies, hard-core potentials give less attrac- 
tive mean- field. However the total self-energy ReS(p, uj) is similar [TTpl . At 
saturation density the dispersive part is not modified much by the three-body 
forces, only the Hartree-Fock energy is lowered. On the other hand at p = 3po, 



Interactions yielding similar binding energies must give similar single-particle en- 
ergies by thermodynamic consistency 
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where the scattering is more important, the dispersive part of the potential is 
lowered when including three-body forces, but the mean-field is similar. De- 
pending on the density, the shift of the single-particle energy caused by the 
three-body forces manifests itself in a different way. In all cases the total self- 
energy Re£(jo, Up) is more attractive when three-body forces are taken into 
account. 

We calculate for the first time the properties of nuclear matter in the self- 
consistent T-matrix approximation with a phenomenological density depen- 
dent three-body term in the nucleon-nucleon interaction. The additional terms 
in the interaction and in the mean-field energy allow to reproduce the empir- 
ical saturation point of symmetric nuclear matter. We find a slightly stiffer 
equation of state at densities above 2p than other approaches (BHF and vari- 
ational). We calculate also the spectral function and the in-medium nuclear 
self-energy. The spectral function shows a stronger quasi-particle peak and 
the potential for the nucleons is more attractive when three-body forces are 
included. The equation of state and the nucleon potential are not directly sen- 
sitive to the details of the short range part of the nuclear interaction. On the 
other hand the proton spectral function at high energies could be a probe of 
short range nuclear correlations. For the interaction studied in this paper we 
find less scattering when density dependent three body terms are included. 
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